Irreversible flow of vortex matter: polycrystal and amorphous phases 
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We investigate the microscopic mechanisms giving rise to plastic depinning and irreversible flow in 
vortex matter. The topology of the vortex array crucially determines the flow response of this system. 
To illustrate this claim, two limiting cases are considered: weak and strong pinning interactions. In 
the first case disorder is strong enough to introduce plastic effects in the vortex lattice. Diffraction 
patterns unveil polycrystalline lattice topology with dislocations and grain boundaries determining 
the electromagnetic response of the system. Filamentary flow is found to arise as a consequence 
of dislocation dynamics. We analize the stability of vortex lattices against the formation of grain 
boundaries, as well as the steady state dynamics for currents approaching the depinning critical 
current from above, when vortex motion is mainly localized at the grain boundaries. On the contrary, 
a dislocation description proves no longer adequate in the second limiting case examined. For 
strong pinning interactions, the vortex array appears completely amorphous and no remnant of the 
Abrikosov lattice order is left. Here we obtain the critical current as a function of impurity density, 
its scaling properties, and characterize the steady state dynamics above depinning. The plastic 
depinning observed in the amorphous phase is tightly connected with the emergence of channel-like 
flow. Our results suggest the possibility of establishing a clear distinction between two topologically 
disordered vortex phases: the vortex polycrystal and the amorphous vortex matter. 



I. INTRODUCTION 

The problem of irreversible [plastic) flow in condensed 
matter transport has attracted a great deal of attention 
over the past decades, due to its striking ubiquity. A di- 
verse variety of physical systems exhibit features of plas- 
tic flow, including vortex matter in type II superconduc- 
tors P, 01, charge density waves 3, supcrfluid ^Hc [4] 
and colloidal matter Q. Plastic transport phenomena 
usually emerge in driven disordered media. The collec- 
tive response of such systems is well described by pin- 
ning theories [1,0, [1| in the limit of weak disorder, where 
dynamics exhibits glassy properties due to the complex 
balance between elasticity of the medium and disorder 
of the substrate. However, as soon as this balance is 
broken, the medium cannot behave elastically and has 
to deform plastically in order to accommodate strain ac- 
cumulation. This is commonly observed whenever the 
strength of disorder overpowers lattice elasticity of the 
medium, or when the medium itself is too soft or di- 
lute to adjust elastically to impurities. In that case it 
may either develop localized topological defects such as 
lattice dislocations or fall into an amorphous phase, not 
reminiscent of the ordered medium, if disorder is strong 
enough. 

Out of the several systems exhibiting this phenomenol- 
ogy, flux line (vortex) lattices in type II superconduc- 
tors represent an ideal framework of investigation, due 
to their remarkable experimental accessibility and their 
ability to be mechanically tuned by changing the applied 
magnetic fleld H The vortex density is proportional 
to H and so is the shear modulus cge of the vortex lat- 
tice in most cases, so that upon changing the value of 
the field one can control the stiffness and the response of 
the lattice. However, as soon as fluctuations (impurities. 



thermal fluctuations etc.) are able to break lattice order, 
mechanical properties of the now-defected vortex lattice 
change dramatically. 

Vortex motion is commonly induced by external cur- 
rent densities J, while collective vortex velocity is pro- 
portional to the experimentally measured voltage V . As 
a consequence, any change in the mechanical response of 
the vortex array is expected to affect the electrodynamic 
response of the underlying superconductor. In this light, 
anomalies in the electromagnetic behavior are usually sig- 
nals of the onset of plasticity. 

The breakdown of the elastic description for a vortex 
lattice can be observed for instance by looking at the 
critical current Jc, which is thus a plastic threshold to 
vortex motion. Experimental studies of the critical cur- 
rent Jc include measurements of Jc as a function of both 
the field H [11, [ll| and the temperature [1| . All systems, 
however, exhibit memory effects governed by an anneal- 
ing process. The system is initially in a disordered state, 
commonly obtained by rapidly cooling the sample from 
a high temperature vortex-liquid phase. This protocol, 
known as field cooling^ ensures that the systems freezes in 
a metastable disordered state. As large enough currents 
are applied, disorder may partially heal and the critical 
currents measured for the healed state are always lower 
than those of the field cooled state. The annealing time 
is found to diverge as the critical current is approached 
from above [ll[ . In fact it was shown that the two states 
may even coexist [l^, [3 , resulting in the coexistence of 
a moving and a pinned phase. It is now widely accepted 
that such a non-annealed fleld-cooled state should corre- 
spond to a disordered phase. This is in agreement with 
the experimental observation of glassy features such as 
power-law distributed relaxation times . 

The so called peak effect is yet another example of the 
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connection between higher critical currents and incipient 
plasticity p^ . A peak in the critical current accompa- 
nies a disordering transition at high temperature or field. 
Peaks in generalized resistivities however appear to be 
ubiquitous whenever a disordering transition is encoun- 
tered, even in rather different systems such as simulated 
polydisperse colloids 

Current- voltage characteristics provide a direct repre- 
sentation of force-velocity relations for the vortex array 
and allow one to disclose the nature of the plastic depin- 
ning transition as soon as Jc is reached. In agreement 
with the observation of coexistence phenomena, one ex- 
pects the transition to be discontinuous. This is a well 
known experimental result [TH and appears very general 
as it has been observed also in driven colloidal matter Q 
as well as in charge density waves [l^ . The discontinuous 
nature is found to increase for increasing disorder. 

As for the topology of the vortex array in the disor- 
dered phase, two possibilities may arise as we discussed 
above, (i) If the lattice is able to retain part of the 
topological order, it will develop dislocations in order to 
relax stresses. Dislocations will arrange into linear pat- 
terns that will act as grain boundaries, resulting in a vor- 
tex polycrystal. Vortex polycrystals have been observed, 
after field- cooling, in various superconducting materials 
such as NbMo[l7|,[i|], NbSes [laMIIII,!!!, BSSCO [H 
and YBCO [2J|- Plastic flow properties of vortex poly- 
crystals are determined by dislocation dynamics and have 
been numerically investigated in several works in the past 
[H, HE [13, On the theoretical side, dislocation argu- 
ments were put forward in order to explain plastic creep 
m Hil and plastic depinning [H [H [M]. (ii) On the 
other hand, if disorder is too strong, or the lattice too 
soft, the dislocation picture may be no longer adequate 
to describe the system. The emergence of completely 
amorphous topo log ies, channel-like flow at relatively low 
currents [ll|, [3l|, [SJl , and smectic dynamics at high drive 
[s^l, summarize the main signatures of the system's be- 
havior in this second scenario. 

Filamentar y v ortex flow is often observed in simula- 
tions [33, m, m, [13, [11]. Depending on simulation pa- 
rameters, it is found to proceed either in individual well- 
defined channels or in several weakly coupled rivers. Such 
observations have also inspired a recent field theoreti- 
cal formulation of the plastic depinning problem, which 
is supposed to provide a reliable picture of irreversible 
transport whenever layered dynamics occur [s^ . 

In this paper we study the mechanisms of plastic de- 
pinning and plastic fiow in vortex matter by means of nu- 
merical simulations. Two limiting cases are considered, 
(i) The first is that of weak pinning interactions. In this 
case disorder is strong enough to produce plastic effects. 
Lattice topology is distorted by the presence of disloca- 
tion arrays, and vortex matter is arranged in a polycrys- 
talline fashion. Critical currents for this systems have al- 
ready been studied in a previous work [28| . Here we focus 
on the stability of vortex lattices against the formation 
of polycrystalline ensembles, the mechanisms of healing. 



and the steady state dynamics for currents approaching 
the depinning threshold from above, when vortex motion 
is mainly localized at the grain boundaries. Filamentary 
flow is found to arise as a consequence of dislocation dy- 
namics, (ii) The second case examined in the paper is 
that of strong pinning interactions. The dislocation de- 
scription proves here inadequate, as the vortex array is 
now amorphous and no remnant of the Abrikosov lattice 
order is left. We study critical currents as a function of 
impurity densities and steady state dynamics above de- 
pinning. We flnd a tight connection between the nature 
of plastic depinning in this regime and the emergence of 
channel-like flow. 

The paper is organized as follows. Section [Til provides 
a description of the numerical methods adopted through- 
out our work. Section [1111 is devoted to the case of weak 
pinning interactions and Section IIVI to the complemen- 
tary results for strong pinning interactions. A discussion 
of the obtained results is found in Section IVl 

II. NUMERICAL METHODS 

We simulate the zero-temperature dynamics of vor- 
tices in a 2D cross-section of linear size L. Such a length 
is fine-tuned in the x and y directions, in order to allow 
an ideal triangular vortex lattice to fit in the simulation 
box. Vortices are subject to external forces produced by 
the applied current density J and interact with Np ran- 
domly distributed pinning centers, which reproduce the 
effect of oxygen vacancies or other impurities in the un- 
derlying material. Unless otherwise noted, simulations 
are performed in a system of size L — 36A with peri- 
odic boundary conditions. From now on, all lengths will 
be expressed in units of A. In particular, we choose a 
value of ^ = 0.2 for the coherence length, corresponding 
to a value of k = 5 for the Ginzburg-Landau parameter. 
Such a value is common in low-Tc superconducting alloys, 
where vortex interactions are well described by Equation 
[5] in the case of thin-film geometries (see below). 

The dynamics of each vortex i at position are de- 
scribed by a set of Ny Langevin equations of motion of 
the form 

Tdrjdt = ^f™(r,-rj) + 5Zf.p(r.-rp + fi(r,), (1) 

j j 

where F is an effective vortex viscosity. The first term on 
the right-hand side describes vortex mutual interactions 
via the long-range force 

f™(r)=AXi (^M^f, (2) 

where A = $g/(87r^A'^), $0 is the quantized flux car- 
ried by the vortices and Ki is a first order modified 
Bessel function [iO] ■ The parameter A sets the force scale 
throghout the simulations. The second contribution in- 
troduces the attraction exerted by each of the Np point 
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defects on vortices. Point defects are randomly located 
at positions (z = 1, ■ • ■ , Np) within the simulation box. 
They exert pinning forces due to a Gaussian potential of 
the form 



V{r - rP) = Va exp 



(3) 



whose amplitude and standard deviation are Vq and ^ 
respectively. If an external current J(r) is applied to the 
sample, it generates a Lorentz-like force acting on the 
vortices 

fi(r) = ^J{r) X z, (4) 

c 

where c is the speed of light. The coupled equations of 
motion ([T]) are numerically solved with an adaptive step- 
size fifth order Runge-Kutta algorithm, imposing peri- 
odic boundary conditions in both directions. 

Transport properties in our simulations are quantified 
by looking at the critical current Jc of the vortex array, 
the collective velocity of the system in the steady state v 
(more precisely, its component along the drift direction) , 
and the probability distribution of velocities in the steady 
state p{v). Topology is taken care of by means of Delau- 
nay triangulations and simulated diffraction patterns. 

Simulations are performed as follows: i) an initial state 
of Ny interacting vortices and Np fixed pinning centers is 
given; ii) we let the system relax until the collective ve- 
locity of the vortex array permanently reaches a value of 
zero; iii) starting form this configuration, different values 
of external driving currents J are tested until the critical 
value Jc is identified. Currents which comply with the 
condition J > identify the region of the phase dia- 
gram where, after an initial transient, the steady state is 
observed with non-zero average vortex velocity. 

We distinguish between two possible simulation proto- 
cols, depending on the initial state that we let relax. In 
the first, the initial state is characterized by a random 
distribution of vortices. This corresponds to a high tem- 
perature state and the subsequent relaxation to a rapid 
quench to zero temperature. Such a procedure mimics 
what is done in field-cooling experiments (FC). In the 
second possible protocol, instead, we choose a perfect 
triangular lattice as the initial state. The experimental 
counterpart of this choice is a zero field- cooling experi- 
ment (ZFC) in which the vortex lattice has been built 
after the temperature quench. 

As mentioned above, our results are organized in two 
main sections. In Section Hill we consider the case of weak 
pinning interactions, choosing a value of Vq = 0.01 for 
the typical depth of the pinning potential wells, as intro- 
duced in Equation [31 In Section |lVl we deal with strong 
pinning interactions, corresponding to Vq = 1; two orders 
of magnitude larger than in the previous case. In both 
cases, the strength of the pinning potential is measured 
in units of the force scale A times the lengthscale A. 
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FIG. 1: Real space images and simulated diffraction patterns 
for the vortex polycrystal. The case of A'^i, = 3120 vortices 
was considered, in interaction with Np = 1032 pinning centers 
(a and c) and Np — 8256 pinning centers (b and d) respec- 
tively. In real space images green, red and blue circles are 
6—, 5— and 7— coordinated vortices respectively. A red-blue 
pair corresponds to a dislocation. A dislocation wall acts as 
a grain boundary, separating grains of different orientations. 
Diffraction patterns were simulated by calculating the struc- 
ture factor in the first Brillouin zone. The peak in the origin 
was removed and lower values were cut off in order to enhance 
clarity. Scale factors are equal in the two plots, so that inten- 
sities in the two realizations can be compared. The patterns 
reproduce multiple rotated replicas of the six-peak structure 
observed for perfect vortex crystals, thus signaling the emer- 
gence of polycrystalline order. For low Np grains are larger (as 
peaks are more intense) and typically less in number, while 
for large Np the grain structure becomes more complex. 



III. WEAK PINNING INTERACTIONS: 
VORTEX POLYCRYSTAL 

In this case, starting form a random vortex configu- 
ration (FC protocol), the system relaxes into a vortex 
polycrystal. The system tries to recover the order of a 
perfect lattice, however impurities counteract this pro- 
cess and the sample freezes into a metastable multi-grain 
configuration, as in Figure [H Lattice order is broken 
by topological defects such as dislocations, which have 
arranged themselves into grain boundaries. When a cur- 
rent above Jc is driven through the sample, the vortex 
polycrystal depins. Depinning is nucleated at the grain 
boundaries and extends to the rest of the lattice as soon 
as large enough currents are applied [28j . 

A detailed study of the topology of the relaxed state 
and its relationship with the critical current is reported in 
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a previous publication ^2§\. For the sake of completeness, 
we recall that the number of impurities Np affects the re- 
laxation process. By varying Np over a range between 64 
and 10320, we observed that all relaxed systems exhibit 
grain structure. However typical grain sizes are found to 
decrease while increasing Np only below a limiting value 
(roughly Np = 4128). Further increases of Np do not 
produce any appreciable drop in grain sizes. We found 
a striking correspondence between this behavior and the 
increase of the critical currents with Np. In particular 
we found that for low values of Np, where grain sizes de- 
crease rapidly, also the critical currents increase rapidly 
(linearly in Np). For large values of Np instead, criti- 
cal currents grow with Np according to a slower square- 
root law. We were able to explain the behavior in terms 
of grain boundary pin ning phenomena. Details can be 
found in Reference |28| . 

Whenever the initial condition is that of a perfect vor- 
tex lattice, as in zero-field cooling experiments (ZFC), the 
system always relaxes into a distorted but not defected 
lattice. The response to a current in such a system is the 
one predicted for elastic depinning [S]. Vortex activity 
is widespread over the sample and it occurs through the 
propagation of avalanches. The critical current is lower 
than in the case of a vortex polycrystal and is propor- 
tional to the density of defects [28| . 



Critical currents and stability of the vortex 
lattice 



Our aim is now to assess the stability of the perfect vor- 
tex lattice state with respect to the formation of a vortex 
polycrystal. To this end, we compared critical current 
characteristics of vortex lattices and vortex polycrystals 
for different vortex densities. The number of vortices in a 
fix-sized sample is proportional to the applied magnetic 
field H. As also the shear modulus of the lattice is pro- 
portional to the magnetic field, one can actually tune the 
mechanical properties of a lattice by changing the vortex 
density. In general one can expect that even for a de- 
fected lattice, the stiffness of the vortex array increases 
as the flux density is increased. Here we show that in- 
creased stiffness also produces increased stability of the 
vortex lattice against the formation of a vortex polycrys- 
tal. 

We start considering the case of high vortex densities 
(iV„ = 3120 vortices in the simulation box), correspond- 
ing to the highest vortex lattice stiffness considered in the 
simulations. We should however emphasize that such a 
density is well below the high-field limit where the dis- 
ordering transition commonly known as peak effect takes 
place. In Figure [5] the critical current as a function of 
the number of defects Np is plotted, both for the vortex 
polycrystal and for the vortex lattice. We observe that 
the two curves remain clearly separated for the whole 
range of pinning densities considered. Indeed, exactly at 
depinning the two states, produced by different sample 
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FIG. 2: Critical current as a function of the number of defects 
for a vortex polycrystal and the corresponding perfect vortex 
lattice in the case of large vortex densities (A'^i, = 3120). The 
two states are well separated and no crossing in found in the 
examined region. Notice that a slightly modified version of 
this figure was first included in a previous publication (Ref. 
[2^ . It is reproduced here for the sake of completeness, as it 
serves a different purpose. 
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FIG. 3; Critical current as a function of the number of de- 
fects for a vortex polycrystal and the corresponding perfect 
vortex lattice in the case of intermediate vortex densities 
{Nv — 2016). The two curves eventually merge meaning 
that the vortex crystal is unstable against the proliferation 
of topological defects and the formation of a polycrystal. 



histories, depin through different phenomena. According 
to the picture given above, motion in the vortex polycrys- 
tal is activated through grain-boundary depinning and 
plastic flow, while the perfect lattice depins elastically 
by releasing widespread avalanches. 

The fact that the two topological states are not equiva- 
lent seems obvious at first sight, but that picture changes 
radically as soon as the vortex density is decreased. In 
Figure [3] we report the same Jc characteristic as above, 
but in the case of intermediate lattice stiffness {Ny = 
2016). The separation between vortex lattice and vortex 
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FIG. 4: Critical current as a function of the number of defects 
for a vortex polycrystal and the corresponding perfect vortex 
lattice in the case of low vortex densities (A*'„ — 1020). The 
two curves almost coincide, except in the case of high purity 
(lowest Np). 



polycrystal is less evident in this case, and for high de- 
fect densities the two states are equivalent at depinning 
in spite of the different sample histories. By looking at 
Delaunay triangulations of the system, one can observe 
that, unlike the case of high vortex densities, here the 
vortex crystal is unstable against the nuclcation of dis- 
locations. The proliferation of topological defects occurs 
for currents just above depinning whenever the impu- 
rity density is high enough. The emergence of a disloca- 
tion network forces the perfect lattice into a state that 
is equivalent to a depinning polycrsytal and this is why 
the two current characteristics eventually merge. This 
phenomenon appears as a re-entrant disordered state in 
the dynamic phase diagram of the vortex lattice for low 
field. Here we provided an interpretation of such a phe- 
nomenon in terms of instability against the formation of 
a vortex polycrystal and showed that for intermediate 
vortex densities it occurs only for large enough impurity 
densities (high Np). 

Finally we consider the case of low vortex densities 
(A^„ = 1020), reported in figure ID This case is emblem- 
atic of lattice instability, as the curve for the vortex sin- 
gle crystal state and the curve for the vortex polycrys- 
tal state merge even at relatively low impurity densities. 
This means the at depinning soft vortex crystals are al- 
most always unstable against the formation of dislocation 
assemblies, unless samples exhibit a high degree of purity. 



B. Nature of the steady state 

So far we have focused on the response of the system at 
depinning, however we have not provided a description of 
the steady state just above depinning. It is known that 
upon applying large enough currents a defected vortex 
lattice might either heal its topological defects and col- 



FIG. 5: Number of topological defects (dislocations) as a func- 
tion of simulation time for high vortex density. The case of 
a vortex polycrystal is considered, with A'^„ = 3120. Two dif- 
ferent curves are shown, for different values of the reduced 
current j = J/Jc — 1. Dotted gray lines are drown as a guide 
to the eye. 



lectively move as a flowing crystal or, if disorder is strong 
enough, or the lattice is particularly soft, it can remain 
topologically disordered and flow in a channel-like fash- 
ion, exhibiting smectic diffraction patterns [35|. These 
considerations apply in the case of high driving currents 
and in a recent numerical study such mechanisms were 
detailed in the case of high AC drive (Ref. [4l|). Here 
instead, we are interested in the behaviour of the system 
for currents close to depinning. 

We can show that in the case of large vortex densi- 
ties {Ny = 3120), a vortex polycrystal always heals into 
an almost perfect vortex lattice, even for currents right 
above depinning. Such a result is reported in Figure [5l 
where the time evolution of the number of defects in a 
single run is recorded for different values of the reduced 
current j = J / Jc— 1. For the largest current considered 
j ~ 1 , the annihilation of dislocations takes place in two 
distinct steps: first a slow process, almost a power-law 
in time, then an abrupt drop in the number of disloca- 
tions, so that only a few pairs survive. Getting closer 
to depinning (j « 0) instead, the process takes a much 
larger time, as experimentally observed in Ref. [llj , and 
goes through a series of steps in the number of defects, 
which suggest that the systems gets momentarily frozen 
in intermediate metastable states. Every time a step is 
overtaken, the collective velocity of the vortex assembly 
experiences a jump, as a decrease in the number of de- 
fects usually implies an increase in the mobility of the 
vortex ensemble. 

In the case of low to intermediate vortex densities, 
however, one expects the above behavior to show sub- 
stantial differences. As a matter of fact, we have already 
shown that even a perfect vortex lattice is in this case 
unstable against the proliferation of dislocations, so one 
can predict that a vortex polycrystal will not heal for 
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FIG. 6: Power spectrum of the collective vortex velocity in 
the drift direction at the steady state just above depinning 
for low vortex densities. In order to obtain a more accurate 
Fast Fourier Transform, we considered the case of = 4032 
vortices and Np = 16512 pinning centers in a simulation box 
of L = 72A, corresponding to the case of low to intermediate 
vortex densities in the usual L — 36A box. Frequencies uj 
are measured in units of (ti — to)"^ ■ No washboard peak is 
observed. One finds a power-law decay instead, reminiscent 
of the 1// noise associated with dislocation dynamics. The 
straight line is drawn for comparison to a law. The 

shape of the cut-off in the power-law behavior is most proba- 
bly due to the nature of the FFT routine. 



currents just above depinning. Indeed, if we plotted the 
time evolution of defects as we did for high vortex den- 
sities, the curve would exhibit large fluctuations but also 
a constant average value in time and no decay, not even 
after very large waiting times. Instead of the time evolu- 
tion of defects, we report the power spectrum of vortex 
velocities in Figure [HI defined as 



k 



(5) 



If the order of a perfect lattice were recovered at least 
to some extent, F(uj) would show a peaked structure 
at the so-called washboard frequency. However in this 
case we observe no peak that is reminiscent of the per- 
fect lattice, meaning that no healing is produced in the 
steady state near depinning for intermediate/low vortex 
densities. In fact, the power spectrum shown in Figure [6] 
has the form of a power-law decaying signal, commonly 
known as broadband voltage noise in experiments [46j. 
This is a well established result, however here we can 
relate it directly to dislocation dynamics. As we know, 
such behavior is due to the fact that the system can- 
not annihilate dislocations (more precisely that it creates 
and annihilates dislocations at the same average rate over 
time). Collective dislocation dynamics is in fact respon- 
sible for the so called 1// noise in plastic deformation of 
crystalline materials, which consists in a power spectrum 
that decays according to a uj~^'^ law [i^, |4^, in good 
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FIG. 7: Steady-state collective vortex velocity in the drift 
direction for A''„ — 1020 vortices as a function of the distance 
to the threshold. Different values of the defect density are 
considered. Inset: logarithmic plot of the same data. 



agreement with our results. As a consequence, the ob- 
servation of broadband voltage noise appears to be again 
a consequence of complex dislocation dynamics in this 
regime. 

We conclude that the phenomenology of the steady 
state near depinning is much richer in the case of low vor- 
tex densities, as no healing is encountered and dislocation 
dynamics heavily affects the response of the system. It 
is known that in this case the vortex drift velocity at the 
steady state approaches zero as the applied current de- 
creases towards Jc in a weakly discontinuous fashion and 
this effect is found to increase in the case of larger pinning 
forces [H, [2^ or alternatively softer vortex arrays [lB| ■ 
In our case however pinning forces are fixed and we are 
able to tune the pinning density instead. The simulated 
steady-state vortex velocities for — 1020 are reported 
in Figure[71 We found that, as expected, the transition is 
indeed weakly discontinuous. The curves collapse as they 
are plotted as a function of the distance to the threshold 
current and moderate variations of the defect concentra- 
tion do not heavily affect the nature of the transition 
in this regime, where dislocation collective dynamics is 
supposed to control vortex flow. Of course, one cannot 
exclude that in principle for higher impurity densities, 
if the vortex ensemble developed amorphous topology, 
a radically different velocity-force relation could be ob- 
served. 

The importance of dislocation dynamics in this regime 
can be also appreciated by looking at the velocity his- 
tograms reported in Figure [8] for the particular case of 
— 1032 vortices and Np = 4128 pinning centers. 
Probability functions for time-averaged individual vor- 
tex velocities in the drift direction at the steady state 
have been calculated for applied currents very close to 
the critical value. Compared to early numerical mea- 
surements of these quantities [1^, here we were able to 
consider larger vortex assemblies and currents closer to 
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FIG. 8: Probability distributions for the components of in- 
dividual vortex velocities in the drift direction. Distribution 
functions are calculated as histograms of time averaged indi- 
vidual velocities in the steady state. Data are represented for 
different values of the reduced current j close to the thresh- 
old {j = 0). Inset: logarithmic plot of the positive tails of the 
above distribution. Velocity distributions are found to decay 
exponentially. 



the threshold. We observe that, very close to depinning 
the probability distribution function has a peak almost 
at zero and a fastly decaying negative tail. The fact that 
a small number of vortices is moving backwards is a direct 
consequence of dislocation motion. If a lattice dislocation 
glides through the array during vortex flow, it produces 
rearrangements of the vortex positions in every spatial 
direction. As a result, the disordered flow of a dislocated 
vortex array in the direction of the applied force requires 
that at every instant of time a small fraction of vortices 
has a negative velocity in the drift direction, if a small 
enough force is applied. This can be observed directly 
by looking at the velocity field under the given condi- 
tions, as reported in Figure M For currents just above 
the threshold, vortex motion is nucleated around moving 
dislocation assemblies and exhibits a whirling dynamics, 
produced by shear stresses associated with grain bound- 
ary glide (Ref. Although the average vortex ve- 
locity has a positive component along the drift direction, 
certain vortices have to accommodate to dislocation dy- 
namics by moving transversely or even backwards. A sim- 
ilar result was obtained in Reference [28*1 , but in the case 
of incipient depinning instead of the steady state. Here 
the low vortex densities ensure that after depinning the 
steady state does not produce global defect annihilation. 
Vortex flow proceeds, mainly along the grain boundaries, 
in meandering filaments. Filamentary flow is thus origi- 
nated by the plastic flow of the vortex polycrystal. 



FIG. 9: Snapshot of the velocity field of a defected vortex 
lattice in the steady state, for applied currents close to the 
critical value. Color scheme is as in Fig. [T] The case of 
low-medium vortex densities is considered, where no healing 
is encountered. Also in the steady state, vortex motion is 
activated around moving dislocation arrays, as in the case 
of incipient depinning (see Reference [IB])- A small fraction 
of vortices accommodates to dislocation glide by moving in 
directions other than that of the applied force (from left to 
right in figure) . 



IV. STRONG PINNING INTERACTIONS: 
AMORPHOUS PHASE 

We now turn our attention to the case of strong pin- 
ning interactions. In Equation [3] we choose Vq = 1 and 
perform the same study proposed in the previous sec- 
tion for the Vo = 0.01 case. The first remarkable feature 
of this regime is that no radical difference is observed 
between the FC and the ZFC protocols. Both a random 
vortex distribution (corresponding to a high temperature 
phase) and a perfect vortex lattice (a zero-field cooled) 
phase relax in the absence of external currents into a 
heavily disordered vortex array, where lattice order is 
completely lost and the dislocation description becomes 
ineffective, as reported in Figure [TOl Such an amorphous 
phase should not be confused with the so-called vortex 
glass phase, induced by high magnetic fields and respon- 
sible for the sudden jump in the critical current known 
as peak effect. In this case, magnetic fields are safely low 
and order loss is due to the strength of impurity fields, 
as commonly observed in heavy-ion irradiated supercon- 
ductors [iH, or even in colloidal aggregates in contact 
with highly heterogeneous substrates Q . However one 
should observe that also in our case we witness an anoma- 
lous current behavior. This aspect will be investigated in 




FIG. 10: Delaunay triangulations,(a) and (b), and simulated 
diffraction patterns, (c) and (d), of relaxed vortex arrays in 
the regime of strong pinning interactions. Data are provided 
for = 3120 vortices in interaction with Np — 2064 point 
impurities, (a) and (c), and Np — 7224 point impurities, 
(b) and (d). In Delaunay triangulations, each circle is an i- 
coordinated vortex, and the following scheme is applied: green 
for i — 6, red for i = 5, magenta for i = 4, blue for i — 7, 
light blue for i = 8 and yellow for other coordination val- 
ues. It appears that the dislocation representation becomes 
inadequate in this the system is found in a purely 

amorphous state. To enhance visibility, vortices are drawn as 
thick circles, however real vortex cores are much smaUer and 
no core overlap effect is produced. Diffraction patterns are 
represented with the same scale factors and cut-offs. Ring 
structures suggest that both translational and orientational 
order are lost. Upon increasing impurity density, the diffrac- 
tion ring becomes thicker and less intense, meaning that the 
vortex spacing distribution has further broadened. 

the following. 

A. Depinning and critical currents 

As in the case of the vortex polycrystal examined 
above, our aim is now to investigate the behavior of the 
amorphous vortex array very close to the critical current 
Jc. As no high-level dislocation theory can be formu- 
lated in this case, one needs to deal with an atomistic 
description of the system at depinning. The simulated 
critical current as a function of the number of impurities 
is reported in Figure [TT] Critical currents are orders of 
magnitude higher than in the case of vortex polycrystals. 
In the range of Ny considered there are no substantial dif- 
ferences between lower and higher vortex densities. The 
three curves start with very low Jc for low impurity den- 
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FIG. 11: Critical current as a function of the number of im- 
purities for the amorphous vortex system in the case of strong 
pinning interactions. Three different values of A*'„ are consid- 
ered. Interpolation lines are drawn as a guide to the eye. 



sities, as if the system was reminiscent of the vortex lat- 
tice, and subsequently experience a sudden jump, where 
the critical current increases severely and reaches a sort 
of plateau. In order to interpret these data correctly, 
we need more information about the nucleation of vortex 
motion at depinning. 

In vortex polycrystals we proved that vortex activ- 
ity was nucleated around depinning grain boundaries (2^. 
The picture that emerges form our simulations is instead 
very different in the amorphous case. Figure [12] shows 
the vortex velocity field at incipient depinning for the 
amorphous vortex array. The system has developed a 
single band of activity, where a few vortices move at 
overly high speed roughly in the direction of the applied 
force (full arrows), surrounded by a stripe of vortices that 
are moving apparently in all possible directions (arrow 
heads alone), while the rest of the sample experiences 
velocities well below the chosen cut-off, thus appearing 
motionless. The interpretation of this picture is rather 
intuitive. In the given regime, pinning forces overpower 
vortex interactions, so that lattice effects are almost neg- 
ligible compared to the case of a polycrystal. However, 
the high forces applied allow the moving vortices to come 
extremely close to the cores of other vortices, experienc- 
ing their strong repulsion. The external force is applied 
uniformly over the sample, however it is able to trigger 
depinning only where it finds a suitable weak spot along 
the force direction. Depinning is then propagated by in- 
dividual vortices, a few units along the given direction. 
Their trajectory is rather heterogeneous, as it is the re- 
sult of the interplay between the high force applied and 
the obstacles encountered in the way: extremely strong 
pinning points and other vortex cores. All around those 
fast vortices, the remaining ones have to rearrange to the 
motion of the former, either being repelled by them or 
trying to fill the empty spots left behind by them. This 
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FIG. 12: Velocity field of the amorphous vortex array at in- 
cipient depinning, for applied currents close to the critical 
value. This time all vortices are represented as green circles 
regardless of their coordination number. Individual velocities 
are represented as arrows. Few vortices have a macroscopic 
velocity and are represented as full arrows. When velocities 
are very small compared to the former, but still above a cer- 
tain cut-off value, they are represented as arrow heads (with 
no line). As velocities fall below the cut-off, no head or line 
is drawn. Depinning proceeds through the fast motion of few 
vortices along a weak direction and the slow wandering of a 
surrounding band. 



corresponds to the band of activity in Figure I12[ which 
becomes broader if a lower velocity cut-off is chosen. 

In a mean-field spirit, each fast vortex moving in the 
drift direction will experience interaction with pinning 
points and surrounding vortices as a sort of effective dy- 
namic friction. This physical picture would then suggest 
that the depinning threshold should be inversely propor- 
tional to the probability p of finding a weak layer in the 
drift direction that would yield vortex motion. Evidently 
p is a function of the number of vortices N^, as in every 
simulation box there are approximately \/Ny layers if dis- 
order is isotropically distributed. Then evidently one has 
that to a first approximation 



Jr oc 



1 



(6) 



This hypothesis is confirmed in FigurefOl where the \/Nv 
scaling of Jc is assessed. Figure [T51 also shows that plot- 
ting the rescaled Jc as a function of Np — Ny, the three 
curves collapse to a master curve in which the current 
jump is located at the origin. The quantity Np — Ny is 
the number of exceeding pinning points or, alternatively, 
the number of free vortices changed in sign. Due to the 
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FIG. 13: Rescaled critical current as a function of the ex- 
ceeding pinning points. Error bars were not represented in 
this plot. The jump in the critical current is observed when 
Np = iV„. 



strength of pinning interactions, vortices are individu- 
ally trapped by impurities. However if Ny exceeds Np, 
at every time step there are Ny — Np vortices that are 
statistically free as they are repelled by locked vortices. 
They are natural weak spots for the system and this ex- 
plains why critical currents are lower in that case. As 
soon as Np reaches Ny, instead, all vortices are locked on 
average, and Jc increases abruptly. Any further increase 
in Np does not affect the system dramatically, thus ex- 
plaining the final plateau region. An estimate for the 
critical current can thus be written as 

Jc = [coQiNy - Np) + cMNp ~ Ny)] (7) 

where cq and ci depend on the pinning potential depth 
Vb and are slowly varying functions oi Np, 9 is the step 
function and the magnetic field dependence is mostly con- 

1/2 

veyed by the Ny factor. To a first approximation one 



has that in this regime 



Jc oc H 



-1/2 



(8) 



B. Steady state and coalescence 



As in the case of the vortex polycrystal examined 
above, we are now interested in the nature of the steady 
state, after depinning is produced. For the range of val- 
ues considered here, we have witnessed no healing or par- 
tial recovery of lattice order. Neither vortex density nor 
larger applied forces seem able to reconstruct an ordered 
lattice. 

If the applied current is very close to the depinning 
value, the steady state does not differ form the picture 
shown in Figure 1121 where few vortices are moving at 
high speed, while the majority are macroscopically inac- 
tive. The coexistence of a moving phase and a pinned 
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FIG. 14: Steady-state collective vortex velocity in the drift 
direction for A'l, — 1020 vortices as a function of the distance 
to the threshold for the amorphous vortex array. Different 
values of the defect density are considered. A logarithmic 
plot is reported in the inset. 



phase can be appreciated in Figure 1141 where velocity- 
force relations are plotted for different impurity concen- 
trations. The universal behavior found in the case of the 
polycrystal is absent here and the strongly discontinuous 
nature of the transition suggests coexistence phenomena, 
as commonly observed in experiments on vortex matter, 
colloids and charge density waves suit- 
able way to quantify coexistence phenomena is to look at 
vortex velocity distributions. In Figure [151 probability 
distribution functions for individual vortex velocities in 
the force direction are represented. As already pointed 
out in a previous numerical study [25j , such distributions 
are bimodal, pointing directly at the two-phase struc- 
ture of the depinning system. However, in our case we 
looked at velocity distributions obtained extremely close 
to the threshold and for a larger system. This allows 
us to access certain quantitative aspects that could not 
be appreciated if such conditions were not met. Figure 
1151 shows two peaks. The first one is around zero veloc- 
ity and quantifies the fact that most vortices are almost 
still. The second peak is found at very high velocities 
and is significantly lower, meaning that only a few vor- 
tices are moving fast. Very close to the threshold (lower 
data) the peak becomes up to 3 orders of magnitude lower 
than the peak in zero. This does not come as a surprise, 
because we have seen that out of thousands of vortices 
only few units contribute to the macroscopic dynamics 
if the force is just above the threshold. As we consider 
lower and lower currents, the tail of the peak in zero 
becomes clearer. The inset of Figure [15] shows that it 
has a power-law behavior over close to two decades in 
velocity, meaning that the slow phase exhibits glassy dy- 
namics. Ideally for currents infinitesimally higher than 
Jc and very large systems, the glassy phase would extend 
to the entire vortex array except an infinitesimal fraction 
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FIG. 15: Probability distributions for the components of in- 
dividual vortex velocities in the drift direction for the amor- 
phous vortex system. Distribution functions are calculated as 
histograms of time averaged individual velocities in the steady 
state. Data are represented for different values of the reduced 
current j close to the threshold (j = 0). Inset: logarithmic 
plot of the positive tails of the above distributions. Only two 
data sets are reproduced here to enhance readability. A power 
law with exponent —3/2 is plotted for comparison purposes. 



of vortex matter. Size effects and the discrete nature of 
vortex matter would make these conditions difficult to 
achieve in experiments. However Figure [T5| shows that 
signals of the glassy phase can also be observed at higher 
currents and in realistic systems. 

Compared to the case of a vortex polycrystal for weak 
pinning interactions, we also observe that the negative 
tail of velocity distributions is now a lot less prominent. 
As no lattice collective rearrangements due to disloca- 
tion glide are present here, vortices can temporarily move 
backwards^ but mainly to accommodate to the motion of 
their fast counterparts. 

As soon as larger currents are applied, new activity 
bands add to the existing one. At sufficiently large cur- 
rent they eventually coalesce and the entire vortex as- 
sembly fiows. However, the strength of impurities is still 
prominent, so that the layered dynamics extends to the 
high drive phase, forcing vortices to move along inde- 
pendent channels. This is the well known smectic phase 
observed even in the case of weak pinning, provided that 
the vortex array is dilute enough. A graphical repre- 
sentation of the smectic phase in our peculiar regime is 
reported in Figure [TBI 



V. CONCLUSIONS 

Topological order of driven vortex arrays in type II 
superconductors appears an essential ingredient to un- 
derstand the flow properties of vortex matter. In most 
common conditions, the Abrikosov vortex lattice order is 
broken by the proliferation of topological defects, such 
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FIG. 16: Evidence of smectic dynamics at very high drive 
in the case of strong pinning interactions, (a) Velocity field; 
(b) simulated diffraction pattern. The typical two-peak struc- 
ture of the diffraction pattern accompanies smectic dynamics. 
Activity bands have coalesced and flow proceeds through de- 
coupled channel motion. 



as dislocations, which add new phenomenology to the in- 
tricate electromagnetic response of this physical system. 
In this paper we emphasize this aspect and examine two 
limiting cases of technological interest: i) the case of weak 
pinning interactions, usually due to intrinsic point disor- 
der in the material, and it) the case of strong pinning 
interactions, which can be extrinsically induced, for in- 
stance, by ion irradiation. In the first case, Delaunay 
triangulations and diffraction patterns unveil a polycrys- 
talline lattice topology, with dislocations assemblies af- 



fecting the current response of the system. On the other 
hand, the topological analysis of the second limiting case 
shows a heavily disordered vortex array, where lattice 
order is completely lost and the dislocation description 
becomes inadequate. 

The depinning of a vortex polycrystal is nucleated at 
the grain boundaries and extends to the rest of the lattice 
as soon as large enough currents are applied (2S | . Vortex 
flow occurs in meandering filaments originated by the 
plastic fiow of the vortex polycrystal. We showed that 
soft vortex lattices are eventually unstable against the 
proliferation of topological defects, and both FC and ZFC 
protocols yield equivalent topological structures, for cur- 
rents just above depinning whenever the impurity den- 
sity is high enough. The average vortex velocity along 
the drift direction in the steady state approaches zero as 
the applied current decreases towards the critical cur- 
rent Jc in a weakly discontinuous fashion. Moreover, 
collective dislocation and grain boundary dynamics in 
this polycrystalline phase give rise to a power spectrum 
that decays according to a cj"^'^ law, in good agreement 
with previous results for plastically deforming crystals 
[43 | . This phenomenology is observed right above depin- 
ning in the steady state attained for low and intermediate 
vortex densities, where no healing is encountered and dis- 
location dynamics heavily affects the response of the sys- 
tem. In this regime, dislocation gliding also induces the 
backwards motion of a small fraction of vortices, which 
can be easily detected in the individual vortex velocity 
distribution. 

The response of the amorphous vortex arrays resulting 
from the action of strong pinning interactions cannot be 
simply interpreted in the language of topological defects 
behavior. Here individual vortex behavior results crucial 
to determine the depinning threshold. The amorphous 
vortex phase appears independently of the cooling pro- 
tocol (FC or ZFC) and critical currents are much higher 
than for their vortex polycrystal counterparts. The mag- 
nitude of Jc is inversely proportional to the probability of 
flnding a layer along the drift direction prompt to vortex 

— 1/2 

motion, and thus scales as Ny . Moreover, it experi- 
ences an abrupt increase as soon as the number of pinning 
sites Np equals the number of vortices , and thus when, 
on average, all vortices may be locked to a pinning point. 
Any further increase in Np does not alter appreciably this 
picture, as confirmed by the final plateau in the critical 
current versus pining density curve. In the steady state, 
the average vortex velocity along the drift direction ap- 
proaches zero as the applied current decreases towards 
the critical current Jc in a strongly discontinuous man- 
ner. Both this quantity and the bimodal distribution 
of individual vortex velocities evidence the coexistence 
of two types of vortex phases right above depinning: a 
glassy phase and a highly mobile phase, which sets up 
fiow bands or channels. 

In conclusion, all the results found for these two lim- 
iting pinning scenarios emphasize the possibility of es- 
tablishing a clear distinction between two dynamical and 
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topologically disordered vortex phases: the vortex poly- 
crystal and the amorphous vortex matter. The behav- 
ior of critical currents and their scaling properties, the 
velocity-force relationships, the velocity histograms and 
the flow geometries corroborate this viewpoint. Recent 
experimental findings on colloidal crystals have already 
illustrated the different nature of these two phases Q. 
New experiments of a similar kind performed in vortex 
matter could help to ascertain the existence of a new 
phase boundary in the already rich dynamic phase dia- 
gram of this fascinating physical system. 
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